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We compute binding energies and root mean square radii for weakly bound systems of N = 4 and 
5 identical bosons. Ground and first excited states of an iV-body system appear below the threshold 
for binding the system with N — 1 particles. Their root mean square radii approach constants in the 
limit of weak binding. Their probability distributions are on average located in non-classical regions 
of space which result in universal structures. Radii decrease with increasing particle number. The 
ground states for more than five particles are probably non-universal whereas excited states may be 
universal. 
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Introduction. The Efimov effect was predicted in [l[ 
as the appearance of a series of intimately related excited 
three-body states when at least two scattering lengths 
are infinitely large. These states can appear at all length 
scales and their properties are independent of the de- 
tails of the potentials. This effect has in recent years 
been studied intensively and extended to a wider group 
of physical phenomena, beginning to be known as Efimov 
physics. In general we define Efimov physics as quantum 
physics where Universality and Scale Invariance apply. 
Universality means independence of the shape of the in- 
terparticle potential. Scale Invariance means indepen- 
dence of the length scale of the system. These conditions 
are rather restrictive but a number of systems are known 
to exist within this window |2|-[7( . The great advantage is 
that one theory is sufficient to explain properties without 
any detailed knowledge of the interactions Q. Further- 
more, properties in different subfields of physics are de- 
scribed as manifestations of the same underlying theory. 

Our physical definition of scale invariance originates 
from the halo physics first realized and discussed in nu- 
clei Q, but quickly observed as applicable also to small 
molecules like the helium trimers [fj. This original defi- 
nition of scale invariance, that the concept applies to any 
length scale is obviously continuous as exemplified by nu- 
clei, atoms and molecules. Often the notion of scale in- 
variance is used in a different mathematical sense where 
the spatial extension of the structures in one given sys- 
tem repeats itself in discrete steps like the factor 22.7 for 
identical particles [l[ . This is a result of the independence 
of potential details and here precisely defining our mean- 
ing with the notion of universality. As the concepts can 
be defined in different ways we will use throughout the 
paper this original physical meaning of scale invariance. 
Together, these two concepts constitute our meaning of 
Efimov Physics which to the best of our knowledge has 
been left undefined in all previous publications. 



The range of validity for such a global theory is only 
well described for two and three particles (3. lil Il0j|. For 
N = 4 two states were found in the zero-range, inherently 
universal, effective field model These states also ap- 
peared as universal in finite-range models in connection 
with each Efimov state [HI]. This is in contrast to [l3| 
where the disentanglement of the scales used to regularize 
the three and four-body zero-range Faddeev-Yakubovsky 
equations gives rise to a dependence of the four-body 
ground state on interaction details. Then a four-body 
scale is needed in analogy to the three-body scale ap- 
pearing independently on top of the two-body proper- 
ties. This apparent discrepancy between Refs. [llL H2I ] 
and 13[ is not yet resolved. 

Recently, three experiments evidenced two four-body 
bound states connected to an Efimov trimer jl4l- ll6| in 
accordance with the theoretical predictions of Ref . [12| ■ 
In two of these experiments were also observed deviations 
from universality [15l . Il6| . Surprisingly, the greatest devi- 
ation were observed for large scattering lengths (— > ±00) 
- exactly at the region where universality should apply 
[l6| . This requires a theoretical explanation where some- 
thing should be added in the universal model. 

Very little is known for five particles with complete so- 
lutions containing all correlations as dictated by the in- 
teraction. With specific assumptions about only s-waves 
and essentially no correlations it was concluded in [13] 
that ground state halos cannot exist for N > 3. These as- 
sumptions are rather extreme and could be wrong or only 
partly correct. However, if halos exist they have universal 
structures as the N = 4 states obtained in llj, [12J . These 
results can only be reconciled by wrong assumptions in 
the halo discussion or by impermissible comparison be- 
tween halo ground states and excited states. 

It was concluded in (3] that Efimov states do not exist 
for N > 3 and furthermore for three par ticles exist only 
for dimensions between 2.3 and 3.8 jl9j. However, by 
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restricting to two-body correlations within the TV-body 
system, a series of (highly) excited TV-body states were 
found with the characteristic Efimov scaling of energies 
and radii 2(| 21 1. Whether they maintain their identity 
and the universal character, when more correlations are 
allowed in the solutions, remains to be seen. 

Two limits to the universality are apparent. The first 
appears for large binding energy where the resulting 
small radii locate the system within the range of the po- 
tentials and sensitivity to details must appear. The less 
strict second limit is for excitation energies above the 
threshold for binding subsystems with fewer particles. 
Structures with such energies are necessarily continuum 
states which may, or more often may not, be classified 
as universal states depending on their structures and the 
final states reached after the decay. 

Even for four particles where universality is found 
111 |l2j , a number of questions are still unanswered. For 
five and more particles the information becomes very 
scarce. A novel study claiming universality for ground 
states of a Van der Waals potential has appeared for par- 
ticle number less than 40 [22]. The critical mass is found 
as a substitute for the critical strength, but the computed 
radii at threshold cannot be reliably extracted. 

The purpose of the present paper is to explore the win- 
dow for Efimov physics. We shall investigate the bound- 
aries for universality preferentially leading to general con- 
clusions applicable to systems of TV-particles. We first 
discuss qualitative features and basic properties, then ex- 
tract numerical results for 4 and 5 particles very close to 
thresholds of binding, and relate to classically allowed 
regions. We only investigate Brunnian (TV-body Bor- 
romean) systems [23| where no subsystem is bound. 

Qualitative considerations. For two particles the in- 
finite scattering length corresponds to a bound state at 
zero energy. Variation of 1/a around zero produces ei- 
ther a bound state of spatial extension a or a continuum 
state corresponding to spatial configurations correlated 
over the radius a. 

For three particles the Efimov effect appears, i.e. for 
the same interaction, a — ±oo (1/a = 0), infinitely 
many three-body bound states emerge with progressively 
smaller binding and correspondingly larger radii [l| . The 
ratios of two and three-body threshold strengths for sev- 
eral potentials were derived in [13, HH Ell ■ These thresh- 
olds for binding one state can be characterized by a value 
of 1/a [HI [H| . Infinitely many bound three-body states 
appear one by one as 1/a is changed from the three- 
body threshold for binding to the threshold for two- 
body binding 1/a = 0. Moving opposite by decreas- 
ing the attraction these states one by one cease to be 
bound. They move into the continuum and continue as 
resonances [26]. For asymmetric systems with a bound 
two-body subsystem the three-body bound state passes 
throu gh t he particle-dimer threshold becoming a virtual 
state [27]. This behavior holds even for particles with 



different masses [281 ] . 

All three-body s-wave states from a certain energy and 
up are universal. However, this is not an a priori obvious 
conclusion but nevertheless true because two effects work 
together, i.e., for 1/a = the system is large for the ex- 
cited Efimov states and for finite 1 /a the binding is weak 
and the radius diverges with binding [5]. Both Efimov 
states and weakly bound states are much larger than the 
range of the interaction. The continuous connection of 
these bound states and resonances is therefore also in the 
universal region. 

The recent results for four particles were that each 
three-body state has two four-body states attached with 
larger binding energy 11, 0]. These four-body states 
are both described as having universal features unam- 
biguously related to the corresponding three-body states 
for interactions of both positive and negative scattering 
lengths. Detailed information of structure, correlations, 
and posssible limits to universality are not available. 

The one-to-one correspondence between the two four- 
body states and one three-body state can perhaps be 
extended such that two weakly bound TV-body states ap- 
pear below the ground state of the (TV — l)-body state. 
This seems to be rather systematic for TV-body Efimov 
states obtained with only two-body correlations 0, [2l| . 
If these TV-body Efimov states remain after extension of 
the Hilbert space to allow all correlations, we can ex- 
pect these sequences to be continued to the thresholds for 
binding by decreasing the attraction. However, ground 
and lowest excited states may be outside the universal 
region but the sequences may still exist. In any case the 
scaling properties are different for the TV-body Efimov 
2C 
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The basic reason for the difficulties in finding detailed 
and general answers is related to the fact that the thresh- 
olds for binding are moving monotonously towards less 
attraction with TV 3, HE]. For TV = 2 weak binding 



and large scattering length is synonymous. Already for 
TV = 3 this connection is broken but the weak binding 
still causes the size to diverge [Bj]. The indications are 
that for TV > 3 the size remains finite even in the limit 
of zero binding. 

Basic properties. We consider a system of TV identical 
bosons each of mass m. They are confined by a harmonic 
trap of frequency ui t corresponding to a length parameter 
bf = h/{mwt). The particles interact pairwise through a 
potential V of short range b <C b t . We shall use the gaus- 
sian shape V = Vq exp(— r 2 /b 2 ). The chosen values of 
Vq, b, and m lead to a two-body scattering length a and 
an effective range R e . The solution to the Schrodinger 
equation is approximately found by the stochastic vari- 
ational method [2i|. The results are energies and root 
mean square radii. 

For two-body systems we know that the n'th radial 
moment only diverge at threshold of binding when the 
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angular momentum I < (n + l)/2, see Q. The equality 
sign implies a logarithmic divergence with binding B 2 in 
contrast to the normal power law B l 2 "Z 2 - 1 / 2 . For the 
mean square radius this implies divergence for I < 3/2. 
For an TV-body system with all contributions entirely 
from s-waves we can generalize these rigorous results 
from two-body systems @. The number of degrees of 
freedom is / = 3 (TV — 1) and the generalized centrifugal 
barrier is obtained with an effective angular momentum 
I* = (f — 3)/2. Divergent root mean square radius is 
then expected when I* < 3/2 or equivalently when / < 6 
or TV < 3. If this result holds, four-body systems should 
have finite root mean square radii even at the threshold 
of binding. 

The size of the system is measured by the square root 
of the mean square radius, ( r 2 /b 2 ), which is expressed in 
units of the "natural" size of the systems, i.e. the range 
of the binding potential. The dimensionless unit of the 
binding energy Bm of the system is B = mb 2 \BN\/h 2 . 
Both Universality and Scale Invariance is therefore de- 
tected by inspection of these quantities as functions of 
parameters and shape of the potentials. In regions where 
the curves are proportional, we conclude that the prop- 
erties are universal and scale invariant. Results for dif- 
ferent potential shapes can be expressed in terms of a 
standard potential by scaling the range. Then the in- 
dividual curves would fall on top of each other in the 
universal regime. 

Clasical allowed region. Universal properties can in- 
tuitively only appear when the structures are outside 
the potentials because otherwise any small modification 
would have an effect on the wavef unction. Consequently 
the property would be dependent on these details in con- 
flict with the assumption of universality. For two-body 
systems the relative wavefunction is therefore universal 
only if the largest probability is found outside the poten- 
tial. This means that this classically forbidden region is 
occupied. The system is extremely quantum mechanical 
and very far from obeying the laws of classical physics. 

To investigate the relation between universality and 
the classical forbidden regions for TV particles we need 
to compare features of universality with occupation of 
classical forbidden regions. For two-body systems this is 
straightforward since the coordinate of the wavefunction 
and the potential is the same. The probability of finding 
the system where the energy is smaller than the potential 
energy is then easy to compute as a simple spatial integral 
over absolute square of the wavefunction. 

For more than two particles the problem is well de- 
fined but the classically forbidden regions (total energy 
is smaller than the potential energy) themselves are diffi- 
cult to locate. We attempt a crude estimate which at best 
can only be valid on average. The energy is computed by 
adding kinetic and potential energy, i.e. 

-B N = N{t 1 } + lN(N-l)(V 12 ) , (1) 



where we choose an arbitrary particle 1 to get the kinetic 
part and a set of particles 1 and 2 to get the potential 
energy. The classical region is defined by having positive 
kinetic energy. For a two-body gaussian potential we 
then obtain an estimate of an average, r c i , for the classical 
radius from 
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If the distance between two particles is larger than r c i 
we should be in the universal region. This value can 
then be compared to the size obtained from the average 
distance between two particles, (r 2 2 ), computed in the 
TV-body system from the mean square radius poj . i.e. 
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Thus in the classical forbidden region r 2 j/6 2 from Eq.® 
should be smaller than (r 2 2 /6 2 ) from Eq.Q. 

The four-body system. We show size versus binding 
energy for TV = 4 in FigJTJ The variation arises by change 
of the strength, Vq, of the attractive gaussian. The sys- 
tem is for numerical convenience confined by an external 
one-body field. However, we are only interested in struc- 
tures independent of that field, i.e. intrinsic properties 
of the four-body system. We therefore increase the trap 
size until the states are converged and located at dis- 
tances much smaller than the confining walls. We now 
know that this happens for four particles in contrast to 
the three-body system where the size diverges when the 
binding energy approaches zero. 

In FigQ] we show results for two trap sizes deviating 
by an order of magnitude and larger than the interaction 
range b = 11.65ao (do is the Bohr radius) by a factor of 20 
and 200, respectively. For large binding in the lower right 
corner the results for the ground state is independent of 
trap size. When the probability extends by about a factor 
of 2 further out than b the effect of the small trap can 
be seen. The tail of the distribution then extends out to 
206 even though the mean square is 10 times smaller. 

In the limit of very small binding energy the radius 
approaches a constant independent of the binding. The 
trap size has to be increased to 2006 before the trap has 
no influence which implies that the probability distribu- 
tion is entirely within that distance when the threshold 
for zero binding is reached. The converged size is about 
2.76 for the ground state. Somewhat surprisingly also 
the first excited state, which also is below the energy 
of the three-body state, has converged to a value, 7.66, 
independent of the trap size. A shape different from a 
gaussian would again lead to constants related through 
specific properties of the potentials, but the ratio would 
remain unchanged. This is precisely as found in two di- 
mensions for three particles 30]. Both states are at the 
threshold on average very much smaller than both traps. 
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FIG. 1: (Color online) The mean square radius as function of 
the four and five-body binding energies, all in dimensionless 
units. The trap sizes for four particles are b t = 230.942ao 
(red lines with 4<), 2630.956ao (black lines with 4>), and 
bt — 372.073ao (blue lines with 5) for five particles. Here 
do is the Bohr radius. We show ground (solid) and excited 
states (long-dashed with particle numbers tagged with an *), 
and "classical" two-body radius (dotted (4) and dot-dashed 
(5)) translated by Eq.©. 



Nevertheless the smallest trap would still influence the 
tail of the distribution. 

In FigfT]we also show the estimated classical average 
distance between pairs of particles within the TV-body 
system. This curve is above the ground state radius for 
large binding. Here the probability is mostly found in 
the classical region within the potential, i.e. in the non- 
universal region. Another potential shape would then 
move these curves. The classical and root mean square 
radius cross each other when the size is slightly larger 
than the range b. This limit for universality is similar to 
the halo condition for universality established in [l7| ■ 
At smaller binding energy the classical radius becomes 
less than the size of the system and the probability is on 
average located outside the potential in the non-classical, 
universal region. 

For the extremely small binding energies close to the 
threshold our estimate of the classical radius diverges log- 
arithmically with binding energy. Thus at some point it 
has to exceed the size of the system which we concluded 
converge to a finite value for zero binding. This is sim- 
ply due to the character of the gaussian potential which 
approaches zero for large radii. Zero energy must then 
be matched by an infinite radius. However, this gaus- 
sian tail is too small to obstruct the convergence of the 
probability distribution to a finite size. This cannot de- 
stroy universality because the tail has no influence on 



the wavefunction in this region far outside the range of 
the potential. For universality only the binding energy 
is decisive as one can see explicitly for the two-body sys- 
tem. For iV-body systems the same result follows from 
the asymptotic large distance behavior of the wavefunc- 
tion expressed in hyperspherical coordinates [3l|. Thus 
the classical average radius argument fails for these ex- 
treme energies when the probability has settled outside 
the range of the short-range potential. 

Five-body system. In FigJT] we also show results for 
N = 5 where convergence is reached for the trap size 
of bt = 372ao- We found two pentamers with energies 
-0.0281-C and -0.0113-C below the four-body thresh- 

old (-0.0103-4-2-) for infinite scattering length. The sizes 
for both ground and excited states increase again with 
decreasing binding energy B§ and approach finite values 
when f?5 = 0. These limiting radii of about 1.946 and 4.0b 
are substantially smaller than corresponding values for 
four particles. Still the largest probability is found out- 
side the potential providing the binding. This strongly 
indicates that also these structures are in the universal 
region. Again their ratio is anticipated to be essentially 
independent of potential shape. This conclusion is sup- 
ported by the comparison in FigQ] to the classical ra- 
dius which always is smaller than the radius of the ex- 
cited state and comparable to the radius of the ground 
state. As argued for four particles the largest binding for 
the ground state corresponds to non-universal structure. 
When the binding energy is about 0.3 in the dimension- 
less units on the figure the universal structure appears. 
This happens at about the same energy as for four parti- 
cles. In both cases the probability is pushed outside the 
potential and universality is expected for smaller binding 
energies. 

Conclusions. We have investigated the behavior of 
Brunnian systems near threshold for binding. Ground 
and first excited state for four and five identical bosons 
appear below the threshold for binding three and four 
particles, respectively. Their radii are for small binding 
energies larger than the range of the potential holding 
them together. The largest part of the probability is 
found in non-classical regions resulting in universal struc- 
tures. For six and more particles the ground states would 
be located inside the potential and thus of non-universal 
structures. Excited states are larger and may still be 
universal but already for seven or eight particles also the 
first excited state is expected to be non-universal. The 
numerical results are obtained for a two-body gaussian 
potential but the features originating from wavefunctions 
in non-classical regions of space are expected to be inde- 
pendent of the potential shape. 
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